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INTRODUCTION 


The  importance  of  the  determination  of  residual  stresses  in  a  prestressed 
thick-walled  cylinder  is  well-known,  and  elastic-plastic  loading  response  has 
been  extensively  studied  (refs  1-5).  Most  of  the  earlier  solutions  for 
residual  stresses  were  based  on  the  assumption  of  elastic  unloading  and  only  a 
few  considered  elastic-plastic  unloading  (refs  2,5).  Bland's  work  (ref  2), 
which  neglects  the  Bauschinger  effect,  is  based  on  the  Treses' s  yield 
criterion  and  the  isotropic  hardening  rule.  Kinematic  hardening  is  the 
simplest  theory  that  can  model  the  Bauschinger  effect  (refs  6,7).  If 
unloading  does  not  occur,  there  is  no  difference  between  the  kinematic  and 
isotropic  hardening  models.  For  unloading  with  reverse  yielding,  the  results 
based  on  these  two  models  will  be  different  as  shown  in  a  recent  paper  (ref  5) 
using  the  ADINA  finite  element  code  (ref  8).  The  von  Mises'  yield  criterion 
and  its  associated  flow  rules  were  used  in  both  models. 

This  report  states  a  closed-form  solution  for  elastic-plastic  loading  and 
unloading  in  pressurized  thick-walled  cylinders  using  Tresca's  yield 
criterion,  its  associated  flow  rule,  and  the  linear  kinematic  hardening  law. 
Numerical  results  are  presented  for  a  closed-end  tube. 

ELASTIC-PLASTIC  LOADING 

Consider  a  thick-walled  cylinder,  internal  radius  a  and  external  radius 
b,  which  Is  subjected  to  internal  pressure  p.  The  material  is  assumed  to  be 
elastic-plastic,  obeying  the  Tresca's  yield  criterion,  the  associated  flow 
theory,  and  a  linear  strain-hardening  rule.  Using  the  isotropic  hardening 

References  are  listed  at  the  end  of  this  report. 


theory,  the  elastic-plastic  solution  has  been  obtained  by  Bland  (ref  2).  In 
order  to  consider  the  Bauschinger  effect,  the  kinematic  hardening  theory  is 
used  here.  Subject  to  the  condition  Oq  >  az  >  Oj.,  Tresca’s  yield  criterion 
for  the  Prager-hardening  rule  (ref  6)  states  that  yielding  occurs  when 

(00-00)  -  (or-o 4.)  -  Kq  (1) 


where 


P  P 

“9  “  c£0  ,  04.  *  cer 

define  the  position  of  the  center  of  the  yield  surface;  c  is  a  material 
constant  and  KQ  the  initial  yield  stress.  The  associated  flow  rule  states 


(2) 


d£0  »  -der  and  de2  *  0  (3) 

For  the  case  of  linear  strain-hardening,  the  yield  stress  curve  can  be 

represented  by  a  straight  line, 

K/K0  *  1  +  neP  and  n  -  (E/Kq)111/ (1-m)  (4) 

P 

where  n  is  a  material  constant  and  the  equivalent  plastic  strain  e  is  defined 
by 

X  /  2 

eP  -  Z2/3  /  {(deQP)2  +  (d erP) 2)  =  —  €0P  (5) 

/T 


The  elastic-plastic  solution  for  the  stresses  and  displacements  can  be 
obtained  explicitly.  The  expressions  in  the  plastic  range  (a  <  r  <  p) 

°r/^o  1  p2  1  p2  p 

=  -  (+1  +  ~)  +  -  nB  (-5  -  1)  -  (l  -  nB)log  - 

ae/K0  2  b2  2  r2 

o2/K0  *  vp2/b2  -  2v(l-nB)  log  -  +  Eez/KQ 

(E/K0)(u/r)  =  (l-2v)(l+v)(or/K0)  -  vEez/Ko  +  (l-v2)p2/r2 


are 

(6) 

(7) 

(8) 


(9) 


and  in  the  elastic  range  (p  <  r  <  b) 


where 


and 


°r/lt°  .  1  (£-  5  fi, 

°e/Ko  2  b2  r2 
oz/K0  -  E£2/Ko  +  vp2/b2 


(E/K0)u/r  -  -  (1+v)  [  p2/r 2  +  (l-2v)p2/b2]  -  vEe*/K, 


z' 


<U-2v) 

E£z/K° “  (b2/!2^!)  (p/Ko) 
p  *  0  (open-end)  ,  1  (closed-end) 


/? 

erl  =  n  +  —  (e/Kq)/ (l- v2) 


(10) 

(11) 

(12) 

(13) 

(14) 


(15) 


The  yield  surface  moves  in  translation  during  plastic  deformation  as  given  by 
a0  “  -ar  *  ( »^3/2 )c£P  and  £P  *  f^P^/r2-!)  (16) 

The  elastic  plastic  surface  p  is  related  to  the  internal  pressure  p  by 

p/K0  -  X-  (l-p2/b2)  +  (l-nE)log  p/a  +  ^  nd  ( p 2/a 2—1 )  (17) 

REVERSE  YIELDING 

If  the  pressure  p  given  by  Eq.  (17)  is  subsequently  removed  completely 
with  no  reverse  yielding,  the  unloading  is  entirely  elastic  and  the  solution 
is  given  by 


°r' 

P  b2 

(18) 

°e' 

i~/iz  - 1 [I  r*  - 11 

(19) 

v  - 

vCor’+ae')  +  Eez' 

(20) 

EV  - 

-(  p-2v)p/(b2/a2-l) 

(21) 

-[(1-v  - 

pv)  +  (l+v)b2/r  2]p/(b  2/a  2-1) 

(22) 

The  residual  stress  system,  which  will  be  denoted  by  two  primes,  is  the  sum  of 


the  system  produced  by  loading  and  that  produced  by  unloading,  i.e.,  or"  =■  or 
+  or' ,  etc.  Assuming  the  kinematic  hardening  rule  and  using  Treses' s 
criterion  subject  to  or"  >  o2"  >  oq",  the  reverse  yielding  will  not  occur  if 

(°r"”ctr)  ~  (O0"-«0)  <  Kq  (23) 

Substituting  the  loading  and  unloading  solution  into  Eq.  (23),  we  can 
determine  the  minimum  pressure  (pm)  for  reverse  yielding  to  occur.  The 
equation  for  Pn  is  given  by 

Pm/Ko  -  (l-a2/b2)  (24) 

Equating  Eqs.  (17)  to  (24),  we  can  determine  the  maximum  amount  of  overstrain 
for  reverse  yielding  not  to  occur. 


ELASTIC-PLASTIC  UNLOADING 


Now  suppose  that  the  loading  has  been  such  that  the  internal  pressure  is 
larger  than  p,,,  given  by  Eq,  (24).  On  unloading,  yielding  will  occur  for  a  <  r 
<  P'  with  p'  <  p.  Using  the  kinematic  hardening  rule  during  unloading  and 
assuming  or"  >  az "  >  a g",  we  have 

(Or’-V)  ~  (O0"-o6-)  -  Kq  (25) 


where 


p 

“r  *  cer  »  o0  *  cee 

Since  the  residual  stress  system  is  the  sum  of  two 


P 

(26) 

systems  produced  by  loading 


and  unloading,  combining  Eqs.  (1)  and  (25)  leads  to 

(Or’-Or')  -  (oe'-oe’)  -  2Kq 


(27) 


(29) 


During  elastic-plastic  unloading,  the  associated  flow  theory  states  that 

P  P  P 

d£g'  »  -der'  <  0  and  dez'  ■  0 

P 

It  has  been  assumed  that  the  sign  of  deg'  is  the  same  throughout  the 

unloading  process  and  that  is  negative.  This  will  be  the  case  when  the 

internal  pressure  is  removed  during  unloading.  Since  cz'  =  cz,e  =  -ez  is 

known,  we  can  use  Hooke's  law  and  the  equilibrium  equation, 

dor*  /  dr  =  (aQ'-ar')/r  (30) 

to  express  oz'  in  terms  of  or' 

az  ~  Eez'  +  2vor'  +  vr(dor'/dr)  (31) 

Since  the  dilation  is  purely  elastic 

(du’/dr)  +  u'/r  +  ez'  =  E“ r(l-2 v)( ar'+a9'+oz' )  (32) 

On  integration  using  Eqs.  (30)  and  (31),  we  obtain 

ru’  =  (l-2v)(l+v)E-1r2  or'  -  vez’r2  +  A  (33) 

where  A  is  a  constant.  The  strain  components  can  be  expressed  in  terras  of  or' 

and  (dor*/dr).  The  plastic  strain  components  are 

P  P 

eg'  *  -er'  =  Ar”2  -  (1-  v2)E“  ^(dOj.1  /dr)  (34) 

Using  Eq.  (31)  together  with  Eqs.  (27)  and  (28),  we  have 

r(dor'/dr)  =  -2(K0-ce9’P)  (35) 

Substituting  Eq.  (35)  into  Eq .  (34)  and  determining  the  constant  A  by  the 
condition  e9'P  *  0  at  p',  we  obtain 

A  =  -2(l-v2)(K0/E)p'  2  (36) 

and 

(E/Ko)e0'P  =  -2(  p'  z/r  2-l )/  [2c/E  +  (l-v2)-!]  (37) 

Equations  (35)  and  (36)  with  the  boundary  condition  at  r  =*  a  suffice  to 
determine  or'  in  the  plastic  region.  The  expressions  for  the  stresses  in  (a  < 
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r  <  p')  are  given  explicitly  by 

or7K0  “  p/K0  -  n8(p' 2/a2-p' 2/r  2)  -  2(l-u8)log(r/a)  (38) 

Oq’/Kq  =  or’/K0  -  2  -  2n8(p,2/r2-l)  (39) 

oz'/K0  “  '»(<’r,+0e,)/Ko  -  Eez/Ko  (4°) 

and  in  (p'  <  r  <  b)  given  by 

0r'/Ko  (41) 

-  ±  (p'2/r2  ?  p’2/b2) 

oe’/K0  (42) 

oz’/K0  =  -2 vp' 2/b 2  -  Eez/K0  (43) 


The  continuity  condition  of  or'  determines  the  relation  between  p'  and  p  as 
given  by 

1  -  p,2/b2  +  2(l-n8)  log  --  +  n«p'  2/a*-l) 

a 

=*  -  (l-p2/b2)  +  (l-u8)log(  p/a)  +  -  nB(p2/a2-l)  (44) 

2  2 

The  yield  surface  moves  in  translation  during  elastic-plastic  unloading 
according  to 

a0'  *  ~“r’  *  CE0,P  (45) 

P 

where  eg'  is  given  by  Eq.  (37). 

NUMERICAL  RESULTS  AND  DISCUSSIONS 

Consider  a  closed-end  thick-walled  cylinder  with  the  following 
parameters:  b/a  =3,  v  =  0.3,  and  E/K0  ■  200.  The  numerical  results  for  the 
displacements,  strains,  and  stresses  during  elastic-plastic  loading  and 
unloading  have  been  calculated.  Figure  1  shows  the  relationship  between  the 
internal  pressure  factor  (p/K0)  and  the  dimensionless  elastic-plastic 
interface  (p/a)  for  various  values  of  the  hardening  parameter,  m  =*  0,  0.05, 
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0.1,  and  0.2.  The  displacements  at  the  inside  and  outside  boundaries  of  the 
tube  (Ua  and  U^)  are  shown  in  Figure  2  as  functions  of  the  elastic-plastic 
interface  for  m  =  0.1.  The  solid  and  dotted  curves  represent  the  displace¬ 
ments  during  loading  and  after  unloading,  respectively.  Figure  3  shows  the 
distribution  of  hoop  stresses  (oq)  during  loading  for  p/a  =  1.0,  1.5,  2.0, 

2.5,  and  3.0.  After  complete  unloading  from  different  stages  of  loading,  the 
corresponding  residual  hoop  stresses  (oq")  are  shown  in  Figure  4.  Reverse 
yielding  occurs  in  a  strain-hardening  tube  with  m  =  0.1  only  when  the  plastic 
portion  (p)  is  larger  than  1.652  a.  In  order  to  show  the  effect  of  hardening 
parameters  (m)  on  the  residual  stress  distribution,  the  numerical  results  are 
presented  in  Figure  5  for  m  =  0,  0.05,  and  0.10.  As  can  be  seen  in  the 
figure,  larger  values  of  hardening  parameter  (m)  tend  to  reduce  the  beneficial 
residual  hoop  stress  at  the  bore. 

All  the  results  presented  in  Figures  1  through  5  are  based  on  the 
kinematic  theory.  We  have  also  calculated  the  results  based  on  the  isotropic 
hardening  theory  (ref  2).  Figure  6  shows  a  comparison  of  two  hardening  rules 
for  the  residual  hoop  stresses  in  a  closed-end  thick-walled  cylinder.  The 
dotted  curves  represent  the  isotropic  hardening  model  with  no  Bauschinger 
effect.  According  to  this  model,  there  is  no  reverse  yielding.  The  solid 
curves  show  the  Bauschinger  effect,  and  reverse  yielding  occurs  in  both  cases, 
p'/a  =  1.098,  1.336,  for  p/a  =  2.0,  3.0,  respectively.  The  residual  hoop 
stresses  at  the  bore  are  0q’7Ko  =»  -0.729,  -0.327  for  p/a  =  2.0,  3.0, 
respectively.  According  to  the  isotropic  hardening  rule,  those  values  of 
0 q"/ko  should  be  -1.060,  -1.318  for  p/a  =  2.0,  3.0,  respectively.  These 
numerical  results  indicate  that  the  effect  of  hardening  rules  on  the  residual 
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hoop  stresses  is  quite  significant,  especially  near  the  bore.  Many  plasticity 
theories  for  reverse  yielding  have  been  proposed  and  reviewed  (ref  9),  and 
many  computer  programs  have  been  developed  (ref  10).  It  is  believed  that  the 
numerical  results  based  on  other  theories  will  fall  within  the  limits  obtained 
by  using  the  kinematic  and  isotropic  hardening  models. 
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Figure  1.  Effect  of  hardening  on  the  relation  between  internal  pressure 
and  elastic-plastic  interface. 
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SMCAR-SCM-0  (PLASTICS  TECH  1 

EVAL  CTR, 

BLDG.  351N) 

SMCAR-TSS  (STINFO)  2 

DOVER,  NJ  07801 

DIRECTOR 

BALLISTICS  RESEARCH  LABORATORY  1 

ATTN:  AMXBR-TSB-S  (STINFO) 

ABERDEEN  PROVING  GROUND,  MD  21005 

MATERIEL  SYSTEMS  ANALYSIS  ACTV 

ATTN :  DRXSY-MP  1 

ABERDEEN  PROVING  GROUND,  MD  21005 


NO.  OF 
COPIES 


COMMANDER 
US  ARMY  AMCCOM 

ATTN:  SMCAR-ESP-L  1 

ROCK  ISLAND,  IL  61299 

COMMANDER 

ROCK  ISLAND  ARSENAL 

ATTN:  SMCRI-ENM  (MAT  SCI  DIV)  1 

ROCK  ISLAND,  IL  61299 

DIRECTOR 

US  ARMY  INDUSTRIAL  BASE  ENG  ACTV 
ATTN:  DRXIB-M  1 

ROCK  ISLAND,  IL  61299 

COMMANDER 

US  ARMY  TANK-AUTMV  R&D  COMD  1 

ATTN:  TECH  LIB  -  DRSTA-TSL 

WARREN,  MI  48090 

COMMANDER 

US  ARMY  TANK-AUTMV  COMD  1 

ATTN:  DRSTA-RC 

WARREN,  MI  48090 

COMMANDER 

US  MILITARY  ACADEMY 

ATTN:  CHMN,  MECH  ENGR  DEPT  1 

WEST  POINT,  NY  10996 

US  ARMY  MISSILE  COMD 
REDSTONE  SCIENTIFIC  INFO  CTR  2 

ATTN:  DOCUMENTS  SECT,  BLDG.  4484 

REDSTONE  ARSENAL,  AL  35898 

COMMANDER 

US  ARMY  FGN  SCIENCE  &  TECH  CTR 
ATTN:  DRXST-SD  l 

220  7TH  STREET,  N.E. 

CHARLOTTESVILLE,  VA  22901 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH  AND  DEVELOPMENT  CENTER, 
US  ARMY  AMCCOM,  ATTN:  BENET  WEAPONS  LABORATORY,  SMCAR-LCB-TL , 
WATERVLIET,  NY  12189,  OF  ANY  ADDRESS  CHANGES. 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST  (CONT’D) 


NO.  OF 
COPIES 


COMMANDER 

US  ARMY  MATERIALS  &  MECHANICS 
RESEARCH  CENTER 
ATTN:  TECH  LIB  -  DRXMR-PL 
WATERTOWN,  MA  01272 

COMMANDER 

US  ARMY  RESEARCH  OFFICE 
ATTN:  CHIEF,  IPO 
P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK,  NC  27709 
COMMANDER 

US  ARMY  HARRY  DIAMOND  LAB 
ATTN:  TECH  LIB 
2800  POWDER  MILL  ROAD 
ADELPHIA,  MD  20783 


DIRECTOR 

US  NAVAL  RESEARCH  LAB 
2  ATTN:  DI'.,  MECH  DIV 

CODE  26-27,  (DOC  LIB) 
WASHINGTON,  D.C.  20375 

COMMANDER 

AIR  FORCE  ARMAMENT  LABORATORY 
1  ATTN:  AFATL/DLJ 

AFATL/DLJG 
EGLIN  AFB ,  FL  32542 

METALS  &  CERAMICS  INFO  CTR 
BATTELLE  COLUMBUS  LAB 
1  505  KING  AVENUE 

COLUMBUS,  OH  43201 


COMMANDER 

NAVAL  SURFACE  WEAPONS  CTR 
ATTN:  TECHNICAL  LIBRARY  1 

CODE  X212 

DAHLGREN,  VA  22448 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH  AND  DEVELOPMENT  CENTER, 
US  ARMY  AMCCOM,  ATTN:  BENET  WEAPONS  LABORATORY,  SMCAR-LCB-TL, 
WATERVLIET,  NY  12189,  OF  ANY  ADDRESS  CHANGES. 


NO.  OF 
COPIES 


1 

1 


1 

1 


1 


